Creating, distributing and freezing entanglement with spin chains 



Irene D'Amico n Brendon W. Lovett and Timothy P. Spiller^tl 

^ Department of Physics, University of York, York YOlO 5DD, United Kingdom 
^ Department of Materials, University of Oxford, 0X1 3PH, United Kingdom 
^ Hewlett-Packard Laboratories, Filton Road, Stoke Gifford, Bristol BS34 8QZ, United Kingdom 

(Dated: February 1, 2008) 

We show how branching spin chains can be used to both generate and distribute entanglement from their 
natural dynamics. Such entanglement provides a useful resource, for example for teleportation or distributed 
quantum processing. Once distributed, this resource can be isolated through mapping or swapping the entan- 
glement into specific qubits at the ends of branches. Alternatively, as we demonstrate for distributed bipartite 
entanglement, applying simple single-qubit operations to the end spin of one or more branches can dynamically 
freeze the entanglement at the branch ends. 
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I. INTRODUCTION 

Over the past few years there has been significant interest in 
the propagation of quantum information through spin chains. 
At a fundamental level, it is interesting to determine just what 
the natural dynamics of these systems permits. However, from 
a more practical quantum technology perspective, the con- 
trolled propagation of quantum information has potential use. 
For macroscopic distances, the current consensus is that pho- 
tons, or other quantum states of light, form the best medium 
for quantum communication, as these can propagate through 
optical fibres or free space with high fidelity'. However, for 
much shorter, microscopic, distances it is possible that other 
media, such as spin chains, could make a very useful contribu- 
tion. For example, in future solid state quantum information 
devices there could be a need to provide quantum commu- 
nication links over microscopic distances, between separate 
quantum processors or registers, or between processors and 
memory, analogous to the conventional communication that 
goes on within the computer chips of today. 

In this paper we use the terminology 'spin chain' to cover 
any physical system whose dynamics can be predicted using 
a Hamiltonian isomorphic to that of a coupled spin chain, 
which we shall define precisely later. In practice, this could 
be strings of actual spins (produced chemically, or fabricated) 
connected through interactions, or it could be a string of quan- 
tum dots or molecules (like fullerenes), containing exciton or 
spin qubits. Strings of trapped atoms form another possibility. 

In order to quantify how well quantum processors or reg- 
isters might communicate, there has been an emphasis on 
how well a quantum state can propagate through a spin 
chain or network. Originally it was shown that a single 
spin qubit could transfer with decent fidelity along a con- 
stant nearest-neighbour exchange-coupled chain^. The fi- 
delity can approach unity if the qubit is encoded into a packet 
of spins'*. Perfect state transfer can also be achieved in more 
complicated systems, with different geometry or chosen un- 
equal couplingsi^i^iiS. Parallel spins chains also enable per- 
fect transfer^''", as does a chain used as a wire, with con- 
trolled couplings at the ends". Related studies have also 
been made on chains of quantum dots'^, chains of quantum 
oscillators'^ '"* and spin chains connected through long range 



magnetic dipole interactionsii. State transfer and operations 
through spin chains using adiabatic dark passage has also been 
proposed"'"'^. 

In this work, we consider a different approach to the prob- 
lem of short distance quantum communication. Rather than 
focus on spin chains providing high fidelity transfer of quan- 
tum information, our approach is to construct spin chain sys- 
tems whose natural dynamics creates high fidelity spatially 
separated entanglement from a simple initial state. We will 
show that this can be achieved in branched spin chain struc- 
tures with certain coupling strengths. The consequences of 
branching were first studied in the context of divided bosonic 
chains, composed of coupled harmonic oscillator&i^ii^ - and 
have since been studied in systems of propagating elec- 
trons^SiSl. In both of these cases, the dynamics of Gaussian 
wave packet type excitations have been investigated. By con- 
trast, our paper deals with the propagation of a single spin-up 
excitation localized on a single site, and how this propagates 
through a network of spin-down states. 

We shall further introduce a technique for: freezing dynam- 
ically generated entanglement. Clearly such entanglement 
is a useful and flexible resource '^-'^-'^'^"'^'-^^'^^'^" ' for quan- 
tum information. For example, it could be used for quantum 
teleportation^^, which enables the transfer of a quantum state, 
or to connect separated quantum registers or processors. One 
advantage of this approach is that decoherence or imperfec- 
tions, which can lead to imperfect entanglement production, 
could be countered by purification techniques^*" prior to use 
of the entangled resource. 

II. SPIN CHAIN FORMALISM 

To introduce our formalism, we first consider a one- 
dimensional chain of N spins, each coupled to their nearest 
neighbours. The Hamiltonian for the system is 
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FIG. 1: A schematic four-site Y spin chain system, with nearest 
neighbour couplings as shown. 



where a\ is the z PauU spin matrix for the spin at site i and 
similarly cr!^ = ct^ ± icr* . For actual spins, Eijl is the local 
magnetic field (in the z-direction) at site i and is the 

local XY coupling strength between neighbouring sites i and 
i + 1. For a coupled chain of quantum dots where each qubit is 
represented by the presence or absence of a ground state exci- 
ton, Ei is the exciton energy and Ji.i+i is the Forster coupling 
between dots i and i + l'^. The spin chain formalism applies 
to both such physical systems, and to any others which can be 
described by an isomorphic Hamiltonian. The computational 
basis notation for the spin states at each site is |0)i = | tz)! 
and |1>, = I U),. 

The total z-component of spin (magnetization), or total ex- 
citon number, is a constant of motion as it commutes with 
H. It is therefore instructive to consider a state of the sys- 
tem consisting of the ground state with the addition of a 
single flipped spin. This state is straightforward to pre- 
pare, assuming local control over a spin at, say, the end of 
a chain, and the flipped spin can be regarding as a (con- 
served) travelling qubit as it moves around under the action 
of the chain dynamics^. An efficient way of representing 
such states in an N spin network is the site basis defined as 
|k) = |0i,02, ...,Ofc_i, lfc,Ofc+i, ...,OAr). A system prepared 
in this subspace remains in it. Now the detailed dynamics de- 
pend on the local magnetic fields or exciton energies, but if 
these are independent of location i then the dynamics favour 
quantum state transfer processes, as already mentioned. We 
also adopt this limit for our work here. 

III. BI-PARTITE ENTANGLEMENT CREATION AND 
DISTRIBUTION 

Rather than using spin chains to transfer quantum states, 
our focus here is to use spin chain systems to prepare spatially 
separated entanglement, which can then be used as a resource 
for teleportation, or providing quantum connections between 
separated quantum processors or registers. The simplest sys- 
tem we consider is a Y structure, chosen to prepare bi-partite 
entanglement from a simple initial state. The smallest defin- 
able Y structure contains four sites, which we label as shown 
in figure [T] 

As seen in the figure, there is a deliberate asymmetry in 
the coupling of the outside sites to the hub of the Y. There is 
one "input" site, coupled to the hub with strength Ji, and two 
"output" sites, coupled to the hub with equal strength J2. The 
Hamiltonian for this system is written: 

Hya = Ji |1) (2| + J2(|2> (3| + |2> (4|) + H.c. (2) 

We proceed by noting that |— ) = 2^^(|3) — |4)) is an eigen- 
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FIG. 2: (a) Example of a qubit network. All the couplings are equal 
and have strength j. (b) Equivalent one-dimensional chain with un- 
equal couplings. 



state of the system. If we initialize our system in the state |1) 
then I — ) is decoupled and plays no part in our system dynam- 
ics. We define the orthogonal state |+) = 2^^ (|3) + |4)) and 
rewrite our Hamiltonian in the {|1) , |2) , |+)} space as: 

i/{.4 = Ji|l>(2| + J2V2|2>(+|+i/.c. (3) 

This shows that our system is equivalent to a ID three-site 
spin chain where the "output" site is the symmetric entangled 
state (10)311)4 + 11)310)4). Such a chain effects perfect state 
transfer from the input site 1 to the entangled output so long 
as the couplings are equal - i.e. if Ji = J2\/2. 

In the following we will show that we can extend this entan- 
glement creation and distribution to spin chain systems with 
longer arms. To do this using the approach introduced above 
is rather cumbersome. Rather, we shall use the method of 
analysis introduced by Christandl et al. in^. They find that 
perfect transfer along a longer chain is possible through the 
use of — 1 unequal couplings along an N site chain, which 
satisfy 

J,,,+i = a^i{N -i). (4) 

where a is a constant that sets the size of the interaction and 
Ji,i+i is the coupling between site i and site i + 1. The authors 
o£^ also show that an array of spins with equal couplings j can 
be projected onto a ID chain of unequal couplings. In partic- 
ular, consider an array that is arranged as a series of columns. 
Let each node of column i be connected to n different nodes 
of adjacent column z + 1, and each node of column z + 1 be 
connected to m different nodes of column i. Christandl et al. 
then show that each column can be projected onto a single 
node of a one dimensional spin chain with coupling jy/mn 
between nodes i and i + 1. An example is shown in figure |2l 
If we analyze the four site Y-chain of figure 1 in this way, it 
is trivial to see that the Y-chain is equivalent to a ID chain 
with couplings Ji, J2\/2, exactly as before. According to 
the general properties of a 3-node chain^, this structure al- 
lows for perfect transfer between the two extremes so long as 
Ji = J2\/2- In this case, the perfect transfer of an excita- 
tion in node 1 is allowed across the Y-shaped structure to the 
sites at the the extremes 3 and 4, where it must be shared be- 
tween the two ends, giving rise to entanglement. Interestingly, 
the same \/2 hub branching factor prevents hub reflection for 
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propagating wave packets in divided chains of harmonic os- 
cillators'^, as well as Gaussian wave packets of propagating 
electrons or magnons22i21. 

We can extend this kind of analysis to larger structures, and 
we now introduce a labelling convention to deal with these. 
Referring to figure [3] we choose to number branches clock- 
wise, starting from the input chain. The length of branch k 
(excluding the hub) is indicated as Ik- A Y-structure is then 
specified by the sequence of numbers {li, I2, Is). Integer la- 
bels i will be used for each of the total N sites in the struc- 
ture as follows: Starting with site i = 1 at the left-most site 
of the input chain, label in sequence up to and including the 
hub. Continue by labelling along output branch 2, to its end, 
and then continue further by labelling branch 3, starting from 
the site adjacent to the hub and going to its end. If the site 
furthest from the hub of branch k is called Uk then ni = 1, 
"2 = ^1 + ^2 + 1 and na = = ?i + /2 + ^3 + 1- 

To demonstrate the generation and perfect transfer of bi- 
partite entanglement to the output spins 712 and 71,3, we have 
performed dynamical simulations using the Hamiltonian of 
([T]i, the condition (|4|i and the branching rule for the couplings 
at the hub spin. Figure |4] shows how the hub branching rules 
applies to the (3, 3, 3) structure, where Ji = Jq — a\/6, J2 = 
J5 = a\/TO, J3 = aVT2, J4 — for perfect state trans- 
fer. Figure |5] shows the result of the simulations for both the 
(3, 3, 3), N = 10, and (10, 10, 10), TV = 31, structures. 
The initial condition is ci = 1 (and all others zero), where Ci 
indicates the amplitude coefficient for the state The up- 
per panel of figure |5] shows the temporal evolution of |cip 
and |cArp for the two systems. We underline that, due to 
symmetry, spins N{— n^) and n2 have the same dynamics 
in this case, and our results show that the excitation is com- 
pletely transferred from site 1 to sites 712 and at time TT/2a 
and periodically returns there at regular time intervals of 7r/a. 
In respect to the rescaled time at, the peak corresponding to 
each revival narrows as the length l^. of the branches increases. 
This is important when considering an optimal practical struc- 
ture for distributing entanglement, since a narrower peak puts 
greater time constraints on any entanglement extraction pro- 
tocol. The lower panel of figure |5] shows the corresponding 
fidelity with respect to the state |+) = 2~2(|n2) + jna)). 

We can quantify the amount of entanglement at the chain 
ends (as a function of time) more formally by using the en- 
tanglement of formation, Ep. This has a value unity for a 
maximally entangled state and zero for an unentangled one. 
Specifically, it measures the number of Bell states required to 




FIG. 5: Results of simulations on a (3,3,3) and a (10,10,10) Y spin 
chain system, with nearest neighbour couplings chosen to satisfy Q 
and the hub branching rule. Initial condition: ci = 1 Upper panel: 
ci|^ and |c]v|^ with respect to the rescaled time at. Lower panel: 
corresponding fidelity with respect to the state |+) = 2^2 (|n2) + 
Ina)). 



create the state of interest and for a two qubit state it is given 
by: 

EAp) = h[^±^), (5) 

wheieh{x) = — a;log2(a::) — (1 — a;) log2(l— a:) is the Shannon 
entropy function, r is the "tangle" or "concurrence" squared: 

r = = [max{Ai - A2 - A3 - A4, 0}]^ . (6) 

The A's are the square roots of the eigenvalues, in decreasing 
order, of the matrix pp ~ p ® p*ay ® , where p* 
denotes the complex conjugation of p in the computational 
basis |00) , |G1) , |10) , lll)^^. 

Ep between qubits n2 ~ 7 and ^ 10 of the (3, 3, 3) 
structure is shown in figure |6l as a function of time following 
initialization of the excitation on site 1 . As would be expected, 
a maximally entangled state is obtained after a time 7r/2a and 
at intervals of ir/a thereafter. 

Maximally entangled states can be produced even if the 
length of the output branches, I2 ~ I3, is different from the 
length of the input branch ^i, as long as the nearest neighbour 
couplings satisfy (|4| and the hub branching rule. As an exam- 
ple in figure|2]we compare the results for the systems (3, 3, 3), 
(5, 2, 2) and (7, 1, 1), with notation and initial condition as in 
figure|5] Our results show the perfect generation of an output 
entangled state, as before, at time t = 7r/2a from a simple 
initial input excitation even for the more asymmetric systems. 
As the asymmetry increases, the width of the peaks decreases, 
but the change is not significant. 
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FIG. 6: Entanglement of Formation of the two chain end qubits 
(n2 = 7 and ns = 10) of the ten site Y chain shown in figure E] 
following the initialization procedure described in the text. 
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FIG. 7: Results of simulations on a (3, 3, 3), (5, 2, 2) and (7, 1, 1) 
Y spin chain system, with nearest neighbour couplings chosen to 
satisfy © and the hub branching rule. Initial condition: ci = 1 
Upper panel: jcip and \cn\'^ as a function of the rescaled time at. 
Lower panel: corresponding fidelity with respect to the state | +) = 
2-4(|n2) + |n3)) 



IV. GENERALIZATION TO p OUTPUTS 

The branching rule we have given for preparation of bi- 
partite entanglement using a Y spin chain structure can be 
extended to the case of the different families of tree-like spin 
chain structures with a single input and n output branches. 
Examples of their members are represented in figure[8] 

Panel (a) corresponds to a 'star-like' family, in which there 
is just one hub and all output branches have the same length. 
Extending the notation early introduced, we can indicate the 
members of this family as (m, I, I, I, I). Following on from 
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FIG. 8: (a) Example of a star-like family member (m, I, I, I, I), (b) 
Example of a bifurcation-like family member. 



the arguments in Section|III] the coupling on the outputs from 
the hub spin are all reduced by a factor of 1/ ^ compared to 
the effective one-dimensional spin chain values. (Again, an 
analogous factor also exists for the condition of no hub reflec- 
tion for a wave packets excitation on the bosonic or propa- 
gating fermion systems discussed previously '^'^°'^'.) The en- 
tangled state which can be created by using this family is a 
W-state, consisting of an excitation shared between p sep- 
arated sites, with the symmetric form 0„3, 0„p) + 
|0„2,1„3, ...,0„J + ... + |0„2,0„3, In real sys- 
tems it is likely to be impractical to prepare p-way entangled 
states using this approach for large values of p. However, 
if three-dimensional physical structures can be built — and 
for example there is the potential for this with quantum dot 
systems — modest values of p > 2, such as 3, 4 and 5 could be 
possible. 

Panel (b) gives an example of what we refer to as a 
bifurcation-like structure. These structures contain more than 
one hub - and for perfect transfer the branching rule must 
be implemented at each. Note that for timing of the dy- 
namical evolution to produce complete excitation transfer to 
the outlying sites, the number of sites along all output paths 
from the initial hub must be the same. The entangled states 
which can be created by using this family share the excita- 
tion between the different separated outlying sites with differ- 
ent unequal weights, e.g. the structure represented in panel 
(b) gives the form |l„2,0„3,0„J/\/2+|0„2,l„3,0„J/2 + 
|0„2 , 0„3 , /2. With only modest branching at each hub, a 
bifurcating structure potentially could be fabricated in a planar 
arrangement. This could in principle give rise to a significant 
number of output branches, whose outlying spins can be en- 
tangled, sharing an excitation in an asymmetric manner, where 
the type of asymmetry depends on the number of hubs and on 
their position. 



V. EXTRACTION OF ENTANGLEMENT 

Having generated and distributed, for example, bi-partite 
entanglement using a Y spin chain system, it is useful to be 
able to isolate this entanglement in some way, so it can be 



5 



used as a resource. Clearly, as illustrated in figure |6] the ar- 
rival time of the entangled state at the output relative to the 
preparation time of the excitation at the input is known. How- 
ever, if the entanglement is being created as a resource (and 
particularly if a number of such states are to be purified be- 
fore use) then it is essential to preserve it once it is created. 
This might be done by extracting it (e.g. with swap opera- 
tions) from the spin chain system and transferring it to a stor- 
age system of qubits, or to members of spatially separated 
quantum registers or processors. See^'* for a method of do- 
ing this with an excitonic quantum dot system. Certainly such 
swapping is likely to be necessary in any excitonic realization 
of a spin chain, if the entangled resource is required to exist 
on timescales longer than the ground state exciton coherence 
times and lifetimes^"*. Another possibility is to physically iso- 
late the two nodes at the end of the chain, after the entangle- 
ment has formed there. If fast local control of the couplings is 
possible, then the state could indeed be frozen at the ends of 
the output branches. To achieve this the couplings have to be 
switched fast on a timescale set by the chain dynamics, as il- 
lustrated by figure|6] and this has to be done simultaneously on 
both output branches. In general both the freezing of the spa- 
tially separated entanglement, and its swapping into longer- 
lived "memory" qubits at the ends of the output branches of 
the spin chain systems, requires fast, coordinated action at the 
ends of both branches. This may be fine if these actions are 
both controlled through external optical pulses. 

However, there is another approach to preserving the entan- 
glement, which could be very effective for some realizations 
of the spin chain systems. This can be implemented when an 
entangled state is (a) an eigenstate of the Hamiltonian ^ and 
(b) it can be reached with just one single qubit local opera- 
tion from the entangled state provided by the natural dynam- 
ics of the system. In this case isolation of the entanglement re- 
quires just a single fast single-qubit operation at one of the two 
branch ends. The simplest example of this is again with the 
four-site system illustrated in figure [T] As already remarked, 
the antisymmetric entangled state |— ) = 2^3 (|3) — |4)) is 
an eigenstate of the system, decoupled from the transfer dy- 
namics (The state of qubits 3 and 4 for a system in the |— ) 
state, after tracing out the others, is the maximally entangled 
2 2(10)311)4— 11)310)4)). After state transfer, our system is 
in the | +) state defined earlier. If a phase flip gate, for example 
|G) |0) and |1) — + — 11), fast on the timescale of the spin 
chain dynamics, is applied to just one chain end qubit (3 or 4) 
of the 1+) state, we obtain |— ). Thus our entangled resource 
can be frozen into an eigenstate of this four-site system. 

Things are not quite so simple in larger Y spin chain sys- 
tems. If a phase flip is applied to the output entangled state 
of a larger Y spin chain system, this opposite parity entan- 
gled state is not an eigenstate of the full system. However, 
the entanglement is effectively "dynamically frozen" by the 
action of the single-qubit phase flip, in the sense that the sub- 
sequent dynamics just involves the spins in I2 and ^3 and the 
|— ) entangled state between spins 712, revives periodically 
and more frequently than if the phase flip is not applied. This 
is illustrated in figure|9]for the systems (3, 3, 3), (5, 2, 2) and 
(7, 1,1), where the results of simulations are presented for the 



initial condition |— ) = 2~^(|n2) — |n3)). Couplings are cho- 
sen to satisfy dU and the hub branching rule. The upper panel 

shows the fidelity with respect to the state |— ) ('Fidelity ^ '). 

These results show the perfect revival of an output entangled 
state with a period shorter than it /a. For the system (3, 3, 3), 
comparison between figure [7] and figure |9] shows that the re- 
vival spacing is 7r/2a in this case. As the dynamics of the 
system is restricted by quantum destructive interference to the 
spins in I2 and I3, the time of revival decreases with decreas- 
ing I2 — h- For the limit case (7, 1, 1) - in which ^2 = ^3 = 1 
- the state |— ) is again an eigenstate of the system, so that the 
entangled state is frozen at the output spins n2 and 713. This 
suggests that a system represented by (m, 1, 1) could be used 
to create entanglement in two qubits on a distant (to spins 
away) array or register In this case the entanglement could 
then be trapped in the two qubits by a single qubit operation 
applied to one of them. This system could also be used as 
a switch, for example to start a computation, or to initialize 
from a distance the creation of entangled states needed for 
measurement-based computation^^. The middle panel shows 
the dynamics of |cip for the spins in ^3 plus the hub spin for 
(3, 3, 3). Clearly the hub is never involved in the dynamics 
(the line labelled as 'a'). The lines labelled from 'b' to 'd' cor- 
respond to the dynamics of spins in Z3 from the hub outward. 
No spin in li is involved in the dynamics (not shown). This 
can be understood from initial condition, in which the dynam- 
ics of the two output branches is out of phase by tt. The sym- 
metry of the structure and of the bonds must therefore gener- 
ate destructive interference in li and the hub. The lower panel 
shows the corresponding dynamics for the (5, 2, 2) system. 
Note that, since in this case I2 ^ h = 2, the dynamics oscil- 
lates perfectly between the spins {n2, n-i] and their respective 
nearest neighbours {712 — 1, 773 — 1}. FigurefTOlshows the de- 
tailed dynamics of the components of c; for {712, 713} (upper 
panel, where the coefficient is purely real) and {712 — 1, 773 — 1} 
(lower panel, where the coefficient is purely imaginary). A 
(maximally) entangled state of the type |— ) is periodically and 
perfectly transferred between the two pairs of spins. For the 
general (a;, 2, 2) system with this initial condition, the dynam- 
ics oscillate between two maximally entangled states, which 
have a phase difference. 

As pointed out above, the destructive interference at the hub 
and in branch li is caused by (i) the antisymmetric initial con- 
dition of the type |— ), (ii) the symmetry of the systems, i.e. 
h = h and (iii) the same sequence of bonds in I2 and ^3. 
What is more interesting is that this is a sufficient condition 
to ensure perfect periodic revival of entanglement at the out- 
put spins. This can be understood by defining an operator S, 
which swaps the two outgoing branches of any Y structure. 
Clearly |+) (|— )) is an eigenstate of S with eigenvalue +1 
(— 1) and, further, S commutes with any Hamiltonian that has 
equivalent outgoing branches. We underline that these need 
not satisfy any special coupling rule, only the condition that 
the two branches have to be the same. S is thus conserved and 
so any state with eigenvalue —1 must always have zero am- 
plitude at the hub and in the input branch. For a finite Y sys- 
tem, an initial state of |— ) will therefore always revive, with 
the hub and input branch having zero amplitude at all times. 
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FIG. 9: Simulations on a (3, 3, 3), (5, 2, 2) and (7, 1, 1) Y spin chain 
system, with nearest neighbour couplings chosen to satisfy l|4j and 
the hub branching rule. Initial condition: |— ) = 2~2(|n2) — Ins)). 
Upper panel: Fidelity with respect to the state |— ) as a function of 
the rescaled time at. Middle panel: \ci\^ for the spins in ^3 plus 
the hub spin for (3, 3, 3). The line labelled as 'a' represents the hub 
dynamics; 'b' to 'd' represent the dynamics of spins in ^3 from the 
closest to the hub (b) to ns (d). Lower panel: jci|^ for spins in ^3 
plus the hub for (5, 2, 2). The line labelled as 'a' represents the hub 
dynamics, 'b' the dynamics of the spin closest to the hub and 'c' the 
dynamics of 713. 



This can be illustrated in numerical examples. In figurefTTIwe 
present the results of simulations on a (3, 3, 3) spin chain, for 
two systems in which the set of nearest neighbour couplings is 
generated using a flat distribution in the interval (0,1], but sat- 
isfying the condition that bonds in I2 and I3 (including the one 
to the hub) have the same sequence. The initial condition is 
given by |— ) = 2^2 (|n2) — ^3)). The upper panel shows the 
fidelity with respect to the state | — ) for two choices of random 
couplings: A (solid line) and B (dashed-dot line) as a function 
of the rescaled time at. As can be seen, the fidelity pattern 
is more complicated than for the case of coupling following 
the prescription given by (|4]i and the hub branching rule. In 
particular peaks are present which do not correspond to maxi- 
mum fidelity. Nevertheless peaks corresponding to fidelity of 
unity appear periodically, with the period determined by the 
particular sequence of bonds. The middle (lower) panel shows 
the dynamics of |cip for spins in I3 and the hub for system A 
(B). We remark that neither system A nor B allows for perfect 
transfer between spin 1 and spins {712, n^}: the condition for 
perfect transfer is a more delicate one and in particular does 
not rely on simple symmetry considerations. 




FIG. 10: Results of simulations on a (5, 2, 2) Y spin chain system, 
with nearest neighbour couplings chosen to satisfy Q and the hub 



branching rule. Initial condition: c„^ 



l/\/2 Dynamics 



of real part of c„2 and c„3 (Upper panel) and imaginary part of c„2 _ 1 
and c„3_i (Lower panel). 



A. Complete Freezing of Dynamically Generated 
Entanglement By Local Operations 

Rather than simply dynamically freezing entanglement by 
shortening its revival time, a combination of the entangle- 
ment extraction ideas can be used to actually freeze spatially 
separated entanglement in systems of the bifurcation fam- 
ily. For example, consider a Y structure where at the end 
of each output branch there is a bifurcation into two final 
spins. This is illustrated in figure [12] In terms of the four end 
spins, the natural dynamics will generate a state of the form 
i(|0, 0,0,1) + |0, 0,1,0) + |0, 1,0,0) + |1, 0,0,0)). This is 
shown in figure [T2]by plotting the branch end spin probabili- 
ties |c„. p as a function of rescaled time at. If at a probability 
peak a phase flip is applied to one spin out of each pair, a state 
oftheformi(|0,0,0, 1)-|0, 0,1, 0) + |0, 1,0,0) -|1, 0,0,0)) 
results. This is an eigenstate of the system and the entangle- 
ment is thus frozen. Although four spins are involved, the spa- 
tial separation achieved is between two pairs of spins, so the 
pairs could each be viewed as (singlet) storage buffers for one 
qubit. In this sense the system contains spatially separated 
and stored bipartite entanglement, which could be released 
for future use by single qubit operations and/or coupling to 
other systems at the branch ends. 



VL CONCLUSIONS 

We have shown how branching spin chain systems can be 
used to both generate and distribute entanglement from their 
natural dynamics. Once distributed, this entanglement can 
be isolated through mapping the entanglement into specific 
qubits at the ends of branches, or, as we have shown for 
distributed bipartite entanglement, applying a simple single- 
qubit operation to one end spin of a branch. Given its simplic- 
ity, this latter approach may be particularly useful in initial 
experimental investigations of such branched systems. 

Distributed entanglement provides a useful resource, for ex- 
ample for teleportation^^ or distributed quantum processing. 
In contrast to the use of spin chains to propagate quantum 
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FIG. 11: Results of simulations of (3,3,3) Y spin chain systems 
with nearest neighbour couplings generated using a flat distribution 
over (0,1], but satisfying the condition that the bonds in I2 and I3 
(and to the hub) have the same sequence. Initial condition |— ) = 
2~2(|n2) — Ina)). Upper panel: Fidelity with respect to the state 
|— ). Middle panel: |cip versus the rescaled time at for spins in Ig 
plus the hub for the random set of couplings A (solid line in upper 
panel). The line labelled as 'a' represents the hub dynamics, and 'b' 
to 'd' the dynamics of spins in h from the closest to the hub (b) to 
ns (d). Lower panel: As for the middle panel but for random set B 
(dashed-dot line in upper panel). Note the different scale on the time 
axis here. 
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FIG. 12: Dynamics of \c„^ \ (for the output branch end spins) ver- 
sus the rescaled time at, for the bifurcation structure shown. This 
follows the initialization condition ci = 1. 



states from one place to another with as high a fidelity as pos- 
sible, there could be some advantage in building up a high 
fidelity entangled resource "offline". Real systems, with their 
inevitable imperfections, will almost certainly degrade trans- 
mission fidelities, even if in principle these approach unity. 
Certainly, with the "offline" resource approach, purification^^ 
could be applied to build up a higher fidelity resource than can 
be achieved by direct transmission. This can then be used to 
transfer quantum states or some form of quantum communi- 
cation. In effect, the concept of a quantum repeater^^ could be 
employed in a solid state, spin chain scenario. 



It is clear from the oscillatory behaviour of the dynamics 
that branched spin chain systems could be operated "in re- 
verse", for example providing a method of detecting certain 
entangled states. If an entangled state is fed into the ends of 
branches it can deliver a signature at a certain single node or 
spin (such as the one we have used as the input node in our 
entanglement creation scenarios). This could be detected with 
just a single-qubit measurement. In more involved scenarios, 
entangled states could also be prepared non-Iocally, through 
a simple conditioning measurement on one or more spins, in 
conjunction with the natural spin chain dynamics acting on 
some simple initial state. 



Finally, we comment that all these related effects result 
from the basic dynamics of the branched spin chain systems 
and simply prepared initial states. Whilst some control is need 
over the couplings to achieve entanglement creation, distribu- 
tion and isolation, there is certainly significant potential for 
branched spin chain systems in solid state quantum process- 
ing and communication. This potential will continue to grow, 
as fabrication or creation of solid state systems that can oper- 
ate as spin chains continues to progress. 
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